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^ [ 1. Introduction 

\ Binary neutron star mergers are expected to be one of the most promising source 

of gravitational waves (GW) for the network of laser interferometric and bar detectors 
becoming operational in the next few years. The merger wave signal is expected to be 
sensitive to the interior structure of the neutron star (NS). The structure of high density 
phases of matter is under current experimental investigation in heavy- ion collisions [ |TJ. 
We investigate the dependence of the merger process and its GW signal on the presence of 
quarks in these phases by performing numerical simulations, where the smoothed particle 
hydrodynamics (SPH) method and the conformally flat approximation for the 3-geometry 
^ . in general relativistic gravity are implemented [ ||. 



2. Equation of state 

Two types of EoS for high density matter, one with quarks and the other without, 
are used in our simulations. In the hadronic model, we describe the stellar matter in 
the high-density region using a EoS based on the relativistic mean field theory with 
the TM1 parameter set using a Lagrangian which contains the low-lying baryons and 
mesons as relativistic fields. The matter at low densities containing a small number of 
protons embedded in the neutron sea with electrons, muons, a-particles and heavy nuclei 
is described using a Thomas- Fermi approximation [||. In the hybrid model, we switch to 
a MIT-bag model based EoS (B=90MeV/fm 3 ) to describe the exotic phases of superdense 
matter with up and down quarks and a negligible small amount of electrons and muons [ 
|J. For the possible mixed phase of hadronic and quark matter which appears for densities 
between 1.8po (po :n uclear saturation density) and 3.7po, we use the multiconserved charge 
phase transition construction [ |l| . In both models, T=0 and /3-equilibrium is assumed 
(see Fig. 0). 
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Figure 1. Pressure-density relation for the 
two EoSs. The two equations coincide up 
to p ~ 5 x 10 14 where the fluid enters the 
quark hadron mixed phase in the hybrid 
case. Above p ~ 10 15 only a pure quark 
phase exists. Note the low adiabatic in- 
dex of the hybrid model EoS in the mixed 
phase of about 1.5. 



Figure 2. Central density- mass relation of a 
series of spherical static NS using both EoSs. 
The maximum mass for the hadronic EoS is 
M ~ 2.21M© while it is M ~ 1.78M in the 
hybrid case. Stars above 'm' have a quark- 
hadron mixed core in the interior, above 'q' 
a pure quark core. 



3. Results 

We start both simulations from irrotational quasi-equilibrium configurations [ [| with 
identical stars each having a gravitational mass in isolation of M = 1.5M Q [ 0. In the case 
of the hybrid model, the component stars have mixed phase matter cores with a radius of 
about 2km. The fluid in the mixed phase has a very low adiabatic index between 1 and 
1.5. Therefore, for constant restmass, a spherical star has a larger central density in the 
hybrid model than in the hadronic model. 

In Fig. ^ we plot the evolution of the maximum density. As expected, the value of the 
hybrid model around the initial phase is slightly larger due to a certain amount of mixed 
phase material. At t ~ 2.5ms, this difference vanishes because the maximum density drops 
below the phase transition density of 1.8po as a result of tidal deformation. Soon after, the 
two stars start merging into a single object. The evolution of the merger objects largely 
depends on the EoS. In the hybrid model, the merger object is contracted further to form 
a high density core because of infalling matter and appearing quark matter which forms a 
compact mixed phase core. Finally the gravitation of the high density core overcome the 
pressure and the centrifugal forces to balance with it and the core starts to collapse. On 
the other hand, the hadronic EoS is much stiffer and the matter is harder to compress. 
Hence, the gravitational contraction is turned into an oscillation of the central object. 

The GW signals (Fig. ||) are almost identical until their amplitudes reach the maximum 
values. Around t ~ 3ms, the amplitude drops down to the minimum as the two component 
stars merge. In the hybrid model, since the merger object contracts faster, the quadrupole 
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Figure 3. The maximum density in units 
of pq. Note the difference during the early 
inspiral phase and the collapse phase. 



Figure 4. Gravitational waveforms of both 
models. The inspiral part is almost identi- 
cal. The collapse of the hybrid model results 
in a faster oscillation. 



moment becomes smaller and the angular velocity becomes larger compared to that of 
hadronic model. This is visible in the wave signal as a smaller amplitude and a phase 
shift. A moment later, when the merger object begins to collapse, the signals begin to 
deviate completely from each other. Due to the higher central density, the hybrid model 
merger object rotates and oscillates faster which results in a higher frequency of the GW 
signal. 

The two models can be distinguished by the GW signal from the merger phase (t > 
2.5ms) whose frequency is above ~ 1kHz. However, in this range, the broadband inter- 
ferometers fail due to shot noise and narrowband resonant bar detectors are needed for 
a successful detection. We are expecting that the merger process of binaries with larger 
masses containing pure quark cores (M > 1.7M , see Fig. 0) will also allow for a distinc- 
tion of the models by the inspiral part of the GW signal which is measurable with the 
broadband interferometer detectors. 
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